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We show that a Calogero-Sutherland type model with anharmonic interactions of 
fourth and sixth orders leads to the matrix model corresponding to the branched polymers. 
We also show that by suitably modifying this model one can also obtain N-particle problems 
which are connected to matrix models corresponding to the pure gravity phase as well as 
corresponding to the transition point between the soap bubble and the branched polymer 
phase. 
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In recent years, Calogero - Sutherland (CS) type N-body problems in one dimension 
have received considerable attention in the literature In particular, remarkable 

connections have been found between such models and seemingly totally different models 
0]. Further, these models are shown to correspond to ideal gas in one dimension with 
fractional exclusion statistics. The purpose of this note is to point out one more such con- 
nection. In particular, we show that the N-body problem with equal mass in 1-dimension 
characterized by (K = 2m = 1, g > — 1/2, B > 0) 

N a2 n 

tf = -£^+£fe^+ B £*. 2 +^I>?> 2 + c (5>?) 3 d) 

i=l 1 i<j 1 i i i 

is related to the matrix model corresponding to the problem of branched polymers. We 
also show that by suitably modifying this model one can also obtain N-body problems 
which are connected to matrix models corresponding to pure gravity phase as well to the 
transition point between the pure gravity and the branched polymer phase. 

Consider a many-body system with the hamiltonian given in eq.(p. The correspond- 
ing Schrodinger equation is given by 

Hip = Eip (2) 

Following Sutherland || we will write the wave function as a product of two wave functions, 
one of which carries the Jastraw factor i.e. the antisymmetric part of the wave function 
and the remaining part which carries the exponential damping terms 

tp = <f>(xi)(p(xi) (3) 

where 

<t> = \{\x i -x j \ x (4) 

i<j 

and 

N N 

^ = exp(-a^x?-/3(^xf) 2 ). (5) 

t=l i=l 

On substituting this wave function in eq.(|2]) we find the relation between the power of the 
Jastraw factor and the strength of the inverse square potential 



A 2 - A = g/2. (6) 
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Solving this for A we get 

A=~[l + (l + 2<7)s]. (7) 

Further, we find that the coefficients in the wavefunction and the coupling constants ap- 
pearing in the hamiltonian are related by 

A = 16a[3, B = 4[a 2 - 0(N + 2 + XN(N - 1))], C = 16/3 2 . (8) 

In case these relations are satisfied then the wave function as given by eq. is the ground 
state eigenfunction and the corresponding eigenvalue is given by 

Eo = -^-t(N + XN(N-1)) (9) 

The ground state eigenfunction of this many body system can be interpreted in terms 
of a matrix model[§]. In particular, square of the modulus of this eigenfunction can be 
interpreted as the weight function of the corresponding matrix integral. For our case it 
reduces to the following matrix integral 

Z N = J dMexp (-^TrM 2 - ^(TrM 2 ) 2 ) (10) 

involving hermitian matrices provided A = 1. Of course, A = 1 means, in the original 
model there in no inverse square interaction. This in turn makes the many body problem 
even simpler. 

This matrix integral has been studied in detail |J . The leading order solution is given 
by taking large iV limit where the eigenvalues scale as 



Xi = VNx(i/VN) = VNx(z) (11) 



with the continuous variable z = i/yN taking values between and 1. In the large N limit 
we can solve the integral in the saddle point approximation. The density of eigenvalues 
u(x) is defined as 

«(x) = % (12) 

and its second moment is 

dxu(x)x 2 . (13) 
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Using these definitions the saddle point equation can be written as 

-x + 2b' cx = P I dy^L. (14) 
2 J-n x-y 

Where P stands for principal value. The eigenvalue range (— //, (jl) is determined by nor- 
malization of the eigenvalue density. The solution to the saddle point equation satisfying 
the proper asymptotics|6[] is 

u{x) = -(I + 2&'c)vV - x 2 . (15) 
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From the normalization as well as the self-consistency conditions flT3| ) it then follows that 

&y = 4 - fi 2 . (i6) 

The free energy of the model in the large N limit now takes the form 



E(b ) = lim IiiZn 



dx-u(x)x 2 + b'c 2 — I iir (l(jtil.i-) a (tj) hi | ./: — // 



(17) 



from where we obtain 



E (b')-E (0) = 1(^ 2 -4) - ^log^ 2 /4 (18) 



It is well known that the matrix models have interpretation in terms of summing over 
random surfaces. The interpretation in our case is that the quartic coupling (3 corresponds 
to the touching random surfaces. This phase, where there are several random surfaces 
touching each other is interpreted geometrically as the branched polymer phase. It is 
well known that in the branched polymer phase of the random surfaces the susceptibility 
exponent is positive and is given by 

7s = 1/2. (19) 

We will not reproduce all the details here but it suffices to say that the analysis of 
can be carried through in a straightforward manner. Following their analysis it is easy 
to see that in our case the susceptibility exponent is indeed 1/2 at b' = —1/16 thereby 
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establishing the connection of our N-body problem with the matrix model for the branched 
polymer phase. 

If we modify the many body potential by adding the term 

N N N N 

v new = a' < + b' 4 + a E A E x ? + D ' E>* - *;) 2 > ( 2 °) 

i=l i=l i=l i=l i<Cj 

and suitably adjust the couplings A' , B', C and D' we can get the susceptibility exponent 
7 S = 1/3 as well as 7 S = —1/2. These cases correspond to the crossover from soap bubble 
phase to the branched polymer phase and the pure gravity phase respectively. For the 
crossover phase, the ground state eigenfunction is again as given by eq. (^) where <f> is as 
given by eq. (f|) while ip is given by 

N N N 

ip = exp(-a J2 x i~ ^E x * ) 2 - 7 E x i) ( 21 ) 

i=l i=l i=l 

where a, (3 and 7 are related to the coefficients in the many body potential as follows 
B = 4[a 2 - 3 7 - 6A 7 (iV - 1) — (3(2 + N + XN(N - 1))] 

A = 16a(3, C = 16/3 2 , A' = 16«7 (22) 

B' = I67 2 , C = 32/3 7 , D' = 4A 7 . 

In this case the ground state energy is as given by eq. (^). 

In the pure gravity case, the ground state eigenfunction is again given by eq. (0) with 
4> being given as before by eq. (|J) while (p is now of the form 

N N 

ip = exp(-a^2x 2 i -7^4) ( 23 ) 

i=l i=l 

where the relation between a, 7 and the many body potential is 

A = C = C = 0, B = 4a 2 - 12 7 (1 -2A(JV-1), 
A' = 16«7, B' = I67, D' = 4A7. 

The ground state energy is now given by 



(24) 



Eo = 2~(d^ {N + N{N ~ 1)X) (25) 

We thus see that the N-body quantum mechanical problem on a line with the hamil- 
tonian ([!]) corresponds to the problem of branched polymers. By adding suitable extra 
potential (p0|), we find that the new many body system corresponds to the pure gravity 
phase as well as the crossover of the soap bubble phase to the branched polymer phase. It 
would be interesting to see if this relation could be further explored to get better under- 
standing of the branched polymer phase. 
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